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This  first  report  of  the  project.  It  incorporates  the  report  on  an  analysis  completed  in  the 
previous  project  on  moving  averages,  variances  and  variograms  for  NIR  from  a  SPOT 
image  of  part  of  Fort  A.  P.  Hill.  In  addition  the  computer  programs  for  these  are  included. 
The  maps  of  the  variances  and  the  diagram  of  the  variograms  show  areas  of  the  image  that 
are  not  stationary  and  where  kriging  is  likely  to  perform  less  well  in  prediction  than  wavelet 
analysis.  The  second  part  of  the  reports  summarises  Dr  Oliver’s  visit  to  the  Topographic 
Engineering  Center  in  June  2000.  The  third  section  comprises  work  on  a  part  of  the  image 
analysed  previously,  called  ‘aphillcut’  to  distinguish  it  from  the  larger  section  of  image  us^. 
This  work  has  explored  in  detail  any  relation  between  elevation  from  a  digital  elevation 
model,  the  raw  elevation  data  and  NIR  and  NDVI.  The  correlations  are  weak  in  spite  of  the 
visual  relation.  The  strongest  relation  was  between  the  moving  averages  for  NDVI  and 
elevation  for  a  window  of  10  x  10  pixels.  The  final  section  describes  the  initial  work  th'St^has 
been  done  on  the  hyperspectral  ‘hymap’  data.  This  is  principally  a  correlation  analysis  and 
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Introduction 


This  is  the  first  report  of  the  current  project.  It  contains  several  sections  of  work.  The  analysis  and 
computer  programs  for  the  first  part  of  the  report  on  the  moving  averages,  variances  and 
variograms  were  completed  in  the  previous  project.  There  was  insufficient  time,  however,  to 
complete  the  written  part  of  this  work.  The  complete  section  is  now  included  although  files  of  the 
results  and  computer  programs  had  been  sent  to  TEC  in  digital  form.  Four  days  of  the  present 
contract  were  spent  working  at  the  Topographic  Engineering  Center  and  one  day  at  the  Virginia 
Institute  of  Marine  Science.  In  addition  time  has  been  spent  on  the  wavelets  and  kriging  paper  for 
submission  to  the  Proceedings  of  the  Sixth  International  Geostatistics  Congress.  This  paper  can 
only  be  accessed  on  CD-Rom  at  present,  the  published  proceedings  will  come  out  in  2001.  The  rest 
of  this  report  will  comprise  work  done  on  the  image  and  digital  elevation  model  for  Fort  A.  P.  Hill, 
and  on  the  latest  hyperspectral  imagery. 


Moving  variance  and  tiled  variograms 

Richard  Webster  and  Margaret  A.  Oliver 

The  computation  of  a  variogram  for  a  region  carries  with  it  the  implication  that  the  underlying 
variation  is  stationary  in  the  intrinsic  sense.  The  analyst  assumes  that  the  expected  differences 
between  places,  at  least  for  small  lag  distances,  are  zero  and  that  the  expected  squared  differences 
are  constant  for  any  given  lag; 

E[Z(x)-Z(x  +  h)]  =  0  (1) 

and  E[{Z(x)-Z(x  +  h)}']  =  2y(h).  (2) 


Here,  in  the  usual  geostatistical  convention,  Z(x)  and  Z(x+h)  are  the  values  of  the  random  variable 
Z  at  positions  x  and  x+h,  where  h  is  the  lag,  and  K**)  is  the  semivariance  at  that  lag. 

The  equations  above  refer  to  the  random  process  Z(x),  and  in  any  one  realization  the  actual  values 
depart  from  its  expectation.  Geostatisticians  are  used  to  this  and  accept  Equation  (1)  without  demur 
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in  many  instances.  Equation  (2)  is  assumed  to  hold  everywhere;  i.e.  the  semivariance  depends  on 
the  lag  only  and  not  on  the  position.  Geostatisticians  are  used  to  this  too;  usually  they  have  too  few 
data  to  explore  its  validity.  Remote  images  from  satellites,  each  comprising  several  hundred 
thousand  pixels,  however,  may  cause  us  to  question  the  assumption  when  some  parts  of  an  image 
are  plainly  more  variable  than  others.  The  SPOT  image  of  AP  Hill  is  one  such;  Figure  1,  a  pixel 
map  of  NIR,  is  an  example.  This  is  the  part  of  A.  P.  Hill  that  we  examined  originally.  With  so  many 
data  we  can  examine  the  changes  in  the  variogram  over  the  region,  and  the  result  can  help  us  to 
decide  whether  our  assumption  of  stationarity  is  reasonable. 

Analysis 

We  computed  three  sets  of  criteria  to  evaluate  the  likely  stationarity;  the  moving  average,  the 
moving  variance,  and  the  local  variogram.  The  first  two  are  defined  as  follows. 

Moving  average.  This  is  computed  by  placing  over  the  image  a  window  containing  n  pixels  and 
summing  the  values  of  NIR  within  the  window: 

^(*c)  =  “Z2(x,).  (3) 

«  i=\ 

where  z(x,)  is  the  value  of  the  ith  pixel  in  the  window,  and  Xc  denotes  the  centre  of  the  pixel.  The 
window  is  moved  one  pixel  at  a  time,  and  at  each  new  position  is  computed.  In  this  way  a 

map  of  moving  averages  is  built. 

Moving  variance.  The  moving  variance,  (o^(Xc),  is  computed  in  an  analogous  way  from 

«-l  /=i 


Again,  a  map  of  the  local  variances  can  be  made. 
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Local  variogram.  Semivariances  can  be  calculated  for  a  window  from 


ri^c) 


1 

2m(h) 


m(h) 


y=i 


(5) 


where  the  index  j  refers  to  those  pixels  for  which  a  paired  comparison  is  possible  at  the  particular 
lag  h,  and  /w(h)  is  the  number  of  such  comparisons. 

Computing  local  variograms  for  the  same  windows  as  for  the  average  and  variance  runs  into  some 
difficulties. 

1.  Experimental  variograms  are  luireliable  with  fewer  than  about  100  data.  So  it  has  seemed 
unreasonable  to  attempt  to  compute  them  for  small  windows,  i.e.  less  than  10  x  10. 

2.  They  are  time-consuming,  though  not  so  time-consuming  as  to  make  the  task  impossible  on  a 
modern  workstation. 

3.  Displaying  the  hundreds  of  thousands  of  variograms  is  fraught. 


So,  instead  of  moving  the  window  one  pixel  at  a  time  we  divided  the  image  into  square  tiles 
without  any  overlap.  Also,  to  provide  an  intelligible  display  we  chose  tiles  of  1 5  x  1 5  —  225  pixels. 
From  experience  we  know  that  this  size  of  sample  estimates  the  local  variogram  well.  We  then 
computed  the  variogram  to  a  maximum  lag  of  12.5  pixels. 

Results 

We  computed  the  moving  averages  and  moving  variances  for  square  windows  of  sides  3,  5,  7,  9, 
11,  and  15  pixels.  The  results  are  displayed  in  Figures  2  to  7.  In  Figure  2  (3  x  3  window)  the 
features  in  the  original  image.  Figure  1  are  still  visible.  For  example  the  lakes  in  the  north  and 
centre  (blue  areas)  and  the  line  of  the  road  running  N-S  (green).  The  latter  is  not  as  evident  as  in 
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Figure  1.  In  Figure  3  (5  x  5  window)  the  lakes  are  still  evident,  but  the  outline  is  no  longer  distinct 
and  the  line  of  the  road  is  barely  visible.  As  the  window  widens  further,  Figures  4  to  7,  the  coarse 
features  of  the  image  become  increasingly  apparent.  The  map  for  the  window  of  9  x  9  pixels  is  very 
similar  to  that  for  the  long-range  structure  from  factorial  kriging  (see  previous  final  report).  There 
is  little  evidence  of  non-stationarity. 

The  margins  of  the  maps,  which  broaden  as  the  window  increases  in  size,  are  the  global  average. 
This  value  is  given  when  the  window  dimensions  cannot  be  fitted  due  to  edge  effects. 

Analogous  maps  of  the  moving  variances  are  shown  in  Figures  8  to  13.  These  tell  a  different  story. 
The  smallest  window,  3x3  =  9  pixels,  shows  local,  mostly  sinuous,  patches  of  image  where  the 
variance  is  much  larger  than  the  general  background.  Figure  8.  In  particular  the  variances  are  large 
at  the  margins  of  the  lakes,  and  their  shape  is  evident.  Increasing  the  window  to  25  pixels  leaves  a 
basically  sinuous  pattern,  but  with  the  patches  of  large  variance  beginning  to  broaden,  Figure  9. 
This  broadening  continues  as  the  window  is  enlarged,  and  the  sinuosity  has  almost  disappeared  and 
is  confined  to  patches  of  medium  variation  with  the  window  of  81  pixels.  The  patches  of  large 
variance  appear  as  blobs,  as  they  do  with  the  wider  windows.  As  for  the  maps  of  the  moving 
averages,  the  uniform  bands  at  the  margins  are  of  the  global  variance. 

Figure  14  shows  the  variograms  in  their  corresponding  tiles,  and  so  is  effectively  a  map.  In  each  tile 
the  experimental  variogram  appears  as  the  set  of  points,  the  left-most  of  which  is  at  the  origin  and 
the  right-most  at  a  lag  of  12.5  pixels  is  scaled  so  that  it  lies  near  the  right  hand  margin  of  the  tile. 
The  ordinates  are  scaled  in  the  range  0  to  20  000  NIR^.  The  horizontal  line  through  the  points  is  the 
local  variance  of  the  data.  The  full  image  is  151  pixels  (columns)  178  (rows).  The  width  of  the  tiles, 
15  pixels,  will  not  divide  exactly  into  either,  and  so  there  are  gaps  at  the  right  and  top 
of  the  image. 

Over  most  of  the  image  the  local  variance  is  small,  and  it  is  difficult  to  see  what  autocorrelation 
there  is.  In  a  relatively  few  squares,  however,  the  variance  is  evidently  strongly  structured.  These 
squares  also  happen  to  be  the  ones  with  the  largest  variances,  and  perhaps  it  is  because  they  have 
the  large  variances  that  it  is  possible  to  display  the  results  in  an  informative  way. 
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Moving  average  for  3  x  3  pixels 
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Figure  2.  Pixel  map  of  the  moving  averages  for  NIR  computed  from  a  square 
window  of  3  X  3  pixels. 
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Moving  average  for  5  x  5  pixels 


Figure  3.  Pixel  map  of  the  moving  averages  for  NIR  computed  from  a  square 
window  of  5  X  5  pixels. 
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Moving  average  for  7  x  7  pixels 
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Figure  4.  Pixel  map  of  the  moving  averages  for  NIR  computed  from  a  square 
window  of  7  X  7  pixels. 
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Figure  5.  Pixel  map  of  the  moving  averages  for  NIR  computed  from  a  square 
window  of  9  X  9  pixels. 
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Figure  6.  Pixel  map  of  the  moving  averages  for  NIR  computed  from  a  square 
window  of  1 1  X  1 1  pixels. 
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Figure  7.  Pixel  map  of  the  moving  averages  for  NIR  computed  from  a  square 
window  of  1 5  x  15  pixels. 
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Moving  variance  for  3  x  3  pixels 
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Figure  8.  Pixel  map  of  the  moving  variances  for  NIR  computed  from  a  square 
window  of  3  X  3  pixels. 
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Moving  variance  for  5  x  5  pixeis 
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Figure  9.  Pixel  map  of  the  moving  variances  for  NIR  computed  from  a  square 
window  of  5  X  5  pixels. 
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Moving  variance  for  7  x  7  pixels 
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Figure  10.  Pixel  map  of  the  moving  variances  for  NIR  computed  from  a  square 
window  of  7  X  7  pixels. 
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Moving  variance  for  9  x  9  pixels 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

T“ 

o 

a> 

00 

h' 

CO 

in 

CO 

04 

LU 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

> 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

CQ 

< 

o 

O) 

CO 

CO 

CO 

04 

Figure  11.  Pixel  map  of  the  moving  variances  for  NIR  computed  from  a  square 
window  of  9  X  9  pixels. 
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Moving  variance  for  11  x  1 1  pixels 


Figure  12.  Pixel  map  of  the  moving  variances  for  NIR  computed  from  a  square 
window  of  1 1  x  11  pixels. 
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Moving  variance  for  15  x  1 5  pixels 


Figure  13.  Pixel  map  of  the  moving  variances  for  NIR  computed  from  a  square 
window  of  1 5  X  1 5  pixels. 
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Conclusions 


The  maps  of  moving  average  show  increasingly  clearly  the  coarse  features  of  the  image  as  the 
moving  window  is  widened.  There  is  little  evidence  of  trend  and  none  sufficient  to  suggest  that  the 
generating  process  is  non-stationary.  The  maps  of  moving  variance,  however,  reveal  small  patches 
where  the  variation  is  much  larger  than  over  most  of  the  image.  They  suggest  that  the  process  is  not 
stationary  locally  in  the  variance  and  that  the  intrinsic  hypothesis  should  not  be  assumed. 

The  steeply  sloping  local  variograms  confirm  this  view.  It  was  precisely  in  these  areas  of  local  non- 
stationarity  that  the  wavelet  analysis  performed  better  than  kriging  in  the  data  reconstruction 
(Oliver  et  al.  2000).  One  thing  that  we  shall  try  later  in  this  project  is  to  use  a  moving  variogram 
for  kriging  to  compare  with  the  wavelet  results. 


Report  on  Visit  by  Dr  Oliver  to  TEC  June  14-21  2000 

Margaret  A.  Oliver 

Part  of  the  first  morning  was  spent  with  Mr  J.  Shine  and  Mr  E.  Bosch  discussing  how  the  time 
should  be  spent  during  the  visit.  Several  things  were  achieved  at  TEC  during  the  four  and  a  half 
days  that  I  was  there. 

1)  Mr  Shine  and  I  started  work  again  on  a  short  paper  for  the  International  Journal  of  Remote 
Sensing.  This  is  now  almost  complete  and  needs  the  final  conclusions  to  be  added.  The  figures 
are  almost  complete;  there  is  still  one  variogram  that  needs  to  be  computed  and  modelled.  This 
paper  will  be  ready  to  be  sent  to  the  Journal  in  October  when  Mr  Shine  visits  England  and  we 
can  finalise  its  contents. 

2)  In  addition  I  went  through  various  Genstat  routines  with  Mr  Shine.  In  trying  to  krige  part  of  an 
image  the  program  gave  us  an  error  message  about  the  grid  size.  Since  my  return  to  England  I 
have  contacted  Rothamsted  Experimrntal  Station  where  Genstat  is  developed  and  I  was  told 
that  this  has  been  resolved  in  the  new  version  of  the  product  which  TEC  is  shortly  to  receive 
when  the  license  is  renewed. 
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3)  In  discussing  the  results  of  comparing  the  digital  elevation  model  with  the  image  data  I 
discovered  that  the  original  data  for  elevation  were  available  without  having  been  smoothed  by 
some  form  of  interpolation.  We  extracted  those  using  Imagine  with  some  help  for  the  area 
where  we  have  been  working.  These  data  have  now  been  re-analysed,  but  there  is  little 
improvement  in  the  correlations. 

4)  Jim  Shine  and  I  were  fortunate  enough  to  have  a  meeting  with  Colonel  Jack  Marin  from  West 
Point.  He  showed  us  some  ways  of  dealing  with  hyperspectral  imagery  using  neural  networks. 
We  have  started  to  analyse  the  ‘hymap’  data  and  we  hope  that  the  variogram  will  identify  bands 
that  describe  different  ground  textures. 

5)  Mr  E.  Bosch  and  I  worked  on  the  soil  data  from  England  and  Wales.  He  has  done  a  wavelet 
analysis  of  a  selected  part  of  the  data  in  the  central  part  of  the  country.  I  did  the  factorial  kriging 
of  the  same  area  to  compare  the  wavelet  analysis  for  the  multispectral  analysis.  We  shall 
continue  with  this  work  when  Mr  Bosch  visits  England  in  August. 

6)  I  spent  one  day  at  the  Virginia  Institute  of  Marine  Science.  Kevin  Slocum  is  registered  as  a  PhD 
student  there  and  I  am  one  of  his  PhD  committee  members.  His  oral  examination  took  place  on 
June  lb***.  It  was  successful  and  from  my  perspective  interesting.  In  the  afternoon  I  gave  a 
lecture  which  included  some  general  geostatistics,  but  mainly  focused  on  the  work  that  I  have 
been  involved  in  at  TEC.  I  presented  the  factorial  kriging  analysis  of  Fort  A.P.Hill  and  also  the 
latest  work  on  wavelets.  The  ensuing  discussion  was  lively. 

One  outcome  of  this  visit  was  a  request  to  teach  a  short  geostatistics  course  there. 
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Fort  A.  P.  Hill:  aphlUcut 


Margaret  A.  Oliver  and  Zoe  L.  Frogbrook 

In  the  previous  final  report  we  examined  the  digital  elevation  model  for  5-m  and  20-m  resolutions. 
Although  there  appears  to  be  a  relation  between  NIR  and  elevation  visually  this  was  not  proved  by 
the  correlation  coefficients  that  we  computed,  which  were  small.  Mr  Shine  and  Dr  Oliver  selected  a 
section  of  the  image  that  we  had  been  working  on  that  minimized  the  areas  of  hard-standing 
because  we  felt  that  these  could  be  reducing  the  overall  correlation. 


Figure  15  a  shows  pixel  map  of  NIR  for  the  part  of  the  image  called  aphillcut:  it  is  in  the  north 
western  part  of  the  original  image,  Figurel.  It  is  an  area  of  76  x  83  pixels  compared  with  the  151x 
178  pixels  of  the  original  image:  it  is  half  the  area.  Figure  15  b  shows  the  pixel  map  of  NDVI 
(Normalized  Vegetation  Index).  The  maps  show  similar  patterns  of  variation.  Variograms  were 
recomputed  and  modelled  for  NIR  and  NDVI,  Figure  16a  and  b,  respectively.  They  show  the 
experimental  semi  variances  as  symbols  and  the  solid  lines  are  the  fitted  models.  Table  1  gives  the 
model  parameters.  They  were  both  fitted  best  by  nested  spherical  functions,  and  their  distance 
parameters  are  very  similar:  the  short-range  component  is  near  to  140  m  and  the  long-range  one  is 
about  230  m.  The  latter  is  shorter  than  that  for  the  larger  image  possibly  because  we  have  fitted  the 
model  to  a  shorter  overall  lag  distance  because  the  experimental  variograms  became  slightly 
irregular. 


Table  1.  Model  parameters  for  selected  properties 


Variable  Model  _ _  Model  parameters _ _ _ 

Co _ _£j _ £2 _ £1 _ £2 _ a 


Sub-sampled  Stable  exponential  0.00 
DEM 

NIR  Double  spherical  0.00 

NDVI  Double  spherical  0.00 


105.7  5.653  1-633 

115.4  109.1  6.577  10.49 

0.002707  0.002187  6.940  11.72 


Note:  Co  is  the  nugget  variance,  c,  and  C2  the  sills  of  the  autocorrelated  variance,  Oi  and  02  the  range  of  spatial 

dependence  and  a  is  the  exponent  parameters  for  the  stable  exponential  model. 


The  original  digital  elevation  model  for  this  area  was  at  a  5-m  resolution.  We  sub-sampled  this  to 
obtain  information  at  the  same  resolution  as  the  SPOT  image  (20-m).  At  the  time  that  these  data 
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Differences 


were  analysed  we  did  not  know  that  the  DEM  comprised  interpolated  data  from  30-m  resolution 
measurements  of  elevation.  Both  the  5-m  and  20-m  data  elevation  data  showed  strong  trend  as  for 
the  larger  area.  The  trend  was  modelled  by  linear,  quadratic  and  cubic  functions  fitted  to  the 
coordinates.  Table  2  shows  the  proportion  of  the  variance  accounted  for  by  the  three  functions  for 
both  data  sets.  Five  transects  along  rows  and  columns  of  the  data  were  selected  from  the  image  for 
a  wavelet  analysis  by  E.  Bosch.  These  were  also  examined  for  trend.  Table  3  gives  the  percentage 
variance  accounted  for  by  trend.  It  is  clear  that  the  trend  is  quite  variable  from  place  to  place. 


Table  2.  Trend  analysis  for  the  original  and  sub-sampled  DEM  data. 


%  variance  accounted  for 

Linear 

Quadratic 

Cubic 

Full  DEM  data 
Sub-sampled  DEM* 

15.7% 

14.5  % 

31.7% 

30.2  % 

41.7% 

40.4  % 

’  DEM  data  was  sub-sampled  to  match  the  co-ordinates  of  the  spectral  values  from  the  SPOT  data. 

Table  3.  Trend  analysis  for  transects  taken  from  the  full  DEM  data. 

Transect _ 

Transect  1:  along  y  axis 
Transect  2:  along  y  axis 
Transect  3:  along  x  axis 
Transect  4:  along  x  axis 
Transect  5:  along  y  axis 

The  variogram  was  computed  for  the  20-m  elevation  data  using  the  residuals  from  the  cubic  trend 
function.  Figure  17  a  shows  the  experimental  variogram.  It  is  upwardly  concave  near  to  the  origin 
and  the  usual  models  that  we  fit  were  not  appropriate.  Since  we  consider  the  Gaussian  model  to  be 
unreliable  we  have  fitted  a  stable  exponential  function  whose  equation  is  given  by 

y(h)  =  c|l  -  exp 

where  c  is  the  sill  variance,  h  is  the  lag  and  a  is  an  exponent  that  must  be  less  than  2.  Figure  17  b 
shows  the  experimental  semivariances  and  the  fitted  stable  exponential  model  for  elevation.  Table  1 
gives  the  model  parameters.  The  approximate  working  range  for  the  exponential  function  is  340  m. 


%  of  variance  accounted  for  by  trend 

5.8 

41.5 

16.5 
21.3 
30.1 
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Lag  distance  /  20  m 


a)  Experimental  variogram  for  the  elevation  data  for  ‘aphillcuf ,  b)  the  fitted  stable 
exponential  model  where  the  exponent  was  1 .6. 


To  use  this  variogram  the  kriging  program  had  to  be  amended,  as  we  have  not  used  it  before  for 
kriging.  Kriging  was  done  on  the  residuals  (for  the  20-m  data)  because  the  presence  of  trend 
violates  the  assumptions  of  geostatistics.  After  kriging  the  trend  was  added  back  to  the  estimates 
using  the  residuals.  Figure  18  shows  the  kriged  estimates  of  elevation  for  the  ‘aphillcut’  area.  A 
visual  inspection  of  the  maps  for  NIR  (Figure  15  a)  and  Figure  18  shows  that  the  valley  systems  in 
the  DEM  correspond  with  the  large  NIR  values,  while  the  interfluves  have  intermediate  NIR 
values. 


For  the  correlation  analysis  the  raw  elevation  data  on  a  30  m  grid  were  used,  Figure  19.  They  were 
compared  with  coincident  pixel  information  for  NIR  and  NDVI.  Figure  20.  Table  4  shows  that 
these  correlations  between  the  raw  data  are  small.  Since  we  know  from  experience  that  local  noise 
in  data  can  sometimes  obscure  relations  between  variables  we  computed  the  correlation  coefficients 
between  the  punctually  and  block  kriged  estimates  and  moving  averages  for  the  DEM,  NIR  and 
NDVI.  Figures  21  and  22  show  the  moving  averages  for  these  variables  for  the  10  by  10  window. 
The  strongest  relation,  albeit  still  weak,  is  between  the  DEM  and  NDVI  for  the  moving  average 
computed  with  a  window  of  10  x  10  pixels.  It  is  evident  that  the  visual  relation  between  the 
observed  pattern  in  the  image  data  and  the  DEM  does  not  appear  to  be  as  strong  statistically  as 
visually. 


Table  4.  Correlation  coefficients  for  DEM,  NIR  and  NDVI. 


Property  1 

Property  2 

Correlation  coefficient 

DEM 

NIR 

-0.172 

DEM 

NDVI 

-0.175 

DEM  residuals 

NIR 

-0.034 

DEM  residuals 

NDVI 

-0.042 

Block  kriged  DEM 

Block  kriged  NIR 

0.235 

Block  kriged  DEM 

Block  kriged  NDVI 

0.219 

DEM  moving  average; 
block  side  of  3 

NIR  moving  average; 
block  side  of  3 

-0.197 

DEM  moving  average: 
block  side  of  3 

NDVI  moving  average: 
block  side  of  3 

-0.203 

DEM  moving  average: 
block  side  of  10 

NIR  moving  average; 
block  side  of  10 

-0.315 

DEM  moving  average: 
block  side  of  10 

NDVI  moving  average: 
block  side  of  10 

-0.330 
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Summary 


The  purpose  of  this  detailed  investigation  between  the  DEM  information  and  the  image  data  is  that 
if  it  had  been  strong  it  could  have  been  used  to  restore  compressed  image  information  using 
cokriging.  The  coregionalization  is  too  weak  to  exploit  in  this  area,  but  it  might  be  worth 
considering  elsewhere. 
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Hyperspectral  imagery 


Margaret  A.  Oliver  and  Zoe  L.  Frogbrook 

The  hyperspectral  imagery  obtained  from  the  hymap  scanner,  a  commerical  h3^erspectral  sensor,  in 
May.  It  took  some  time  to  extract  from  the  CD  as  it  was  in  a  format  we  were  not  familiar  with. 
There  are  126  bands  that  have  been  separated  from  the  spectrum.  The  pixel  size  is  7  m  x  7m.  Our 
aim  in  the  analysis  is  to  determine  whether  the  texture  expressed  by  different  bands  can  be 
identified  using  the  variogram,  kriging  and  factorial  kriging  if  the  variograms  are  nested. 

Since  the  data  set  is  so  large  we  sampled  it  and  retained  one  pixel  in  every  three.  This  resulted  in  a 
21-m  separation  between  the  pixel  centres.  Our  preliminary  variogram  analysis  of  the  first  8  bands 
showed  that  these  all  have  a  similar  spatial  structure  to  each  other. 

Table  5  gives  the  correlation  coefficients  between  all  of  the  bands  and  it  is  evident  that  the  first  36 
are  strongly  correlated.  It  is  likely  that  these  will  contain  similar  spatial  information.  The  remaining 
bands  are  moderately  to  weakly  correlated  with  the  first  36  bands,  but  other  groups  of  wavebands 
show  strong  correlation  with  each  other.  Do  the  strongly  correlated  bands  imply  that  spatially  these 
are  mainly  redundant  and  that  we  can  look  at  one  to  obtain  the  information  of  value?  We  shall 
examine  some  of  the  strongly  correlated  bands  to  assess  whether  this  is  so.  We  shall  use  the 
correlation  matrix  as  an  aid  to  selecting  bands  likely  to  contain  information  that  is  different  from 
others. 

Tables  6  and  7give  the  results  of  a  principal  component  analysis  (PCA).  This  multivariate  analysis 
is  one  way  of  reducing  the  amount  of  data  that  is  examined.  Here  we  include  only  the  preliminary 
analysis.  The  first  two  latent  roots  account  for  almost  80%  of  the  variation.  Variograms  were 
computed  from  the  scores  of  the  first  five  components,  that  for  PCS  was  pure  nugget.  This  root 
expressed  only  1.45%  of  the  variation  and  contains  little  or  no  spatial  information.  Nested  models 
provided  the  best  fit  to  the  variograms  of  the  first  two  components  and  exponential  functions  fitted 
PCs  3  and  4  best.  The  model  parameters  are  given  in  Table  7.  The  spatial  scales  for  PCs  3  and  4 
are  similar;  it  is  the  sill  variances  that  are  different.  For  PCs  1  and  2  the  latter  appears  to  have  larger 
ranges  of  spatial  dependence  for  both  the  long-  and  short-range  components. 
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Table  5  (cont.). 


b37 

b38 

b39 

b40 

b41 

b42 

b43 

b44 

b45 

b46 

b47 

b48 

b37 

1.00 

b38 

1.00 

1.00 

b39 

1,00 

1.00 

1.00 

b40 

1.00 

1.00 

1.00 

1.00 

b41 

1.00 

1.00 

1.00 

1.00 

1.00 

b42 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

b43 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

b44 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

b45 

0.99 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

b46 

0.99 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

b47 

1.00 

1.00 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

1.00 

b48 

0.99 

0.99 

0.98 

0.98 

0.98 

0.98 

0.98 

0.98 

0.98 

0.98 

1.00 

1.00 

b49 

0.98 

0.98 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.99 

1.00 

b50 

0.98 

0.98 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0,99 

1.00 

b51 

0.97 

0.97 

0.97 

0.97 

0.96 

0.96 

0.96 

0.97 

0.97 

0.97 

0.99 

0.99 

b52 

0.97 

0.97 

0.97 

0.96 

0.96 

0.96 

0.96 

0.96 

0.97 

0.97 

0.99 

1.00 

b53 

0.97 

0.97 

0.97 

0.97 

0.97 

0.96 

0.96 

0.97 

0.97 

0.97 

0.99 

1.00 

b54 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.98 

0.99 

1.00 

b55 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.98 

0.99 

0.99 

b56 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.98 

0.99 

1.00 

b57 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.98 

0.99 

1.00 

b58 

0.97 

0.97 

0.97 

0.96 

0.96 

0.96 

0.96 

0.96 

0.97 

0.97 

0.98 

0.99 

b59 

0.96 

0.96 

0.96 

0.95 

0.95 

0.95 

0.95 

0.95 

0.96 

0.97 

0.98 

0.99 

b60 

0.94 

0.95 

0.94 

0.94 

0.94 

0.94 

0.94 

0.94 

0.95 

0.95 

0.97 

0.98 

b61 

0.92 

0.93 

0.93 

0.92 

0.92 

0.91 

0.91 

0.92 

0.93 

0.94 

0.95 

0.97 

b62 

0.90 

0.91 

0.91 

0.90 

0.89 

0.89 

0.89 

0.90 

0.91 

0.92 

0.93 

0.94 

b63 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.05 

0.06 

0.06 

0.07 

0.07 

0.07 

b64 

0.11 

0.11 

0.11 

0.11 

0.10 

0.10 

0.10 

0.10 

0.11 

0.12 

0.13 

0.15 

b65 

0.23 

0.24 

0.23 

0.22 

0.22 

0.21 

0.21 

0.22 

0.23 

0.25 

0.28 

0.32 

b66 

-0.12 

-0.13 

-0.13 

-0.13 

-0.13 

-0.13 

-0.13 

-0.13 

-0.13 

-0.13 

-0.12 

-0.12 

b67 

-0.12 

-0.12 

-0.12 

-0.12 

-0.12 

-0.12 

-0.13 

-0.13 

-0.13 

-0.12 

-0.12 

-0.11 

b68 

-0.10 

-0.10 

-0.10 

-0.10 

-0.10 

-0.10 

-0.11 

-0.11 

-0.10 

-0.10 

-0.09 

-0.08 

b69 

0.41 

0.41 

0.41 

0.39 

0.38 

0.38 

0.38 

0.39 

0.40 

0.43 

0.47 

0.53 

b70 

0.44 

0.44 

0.44 

0.42 

0.42 

0,41 

0.41 

0.42 

0.43 

0.46 

0.50 

0.56 

b71 

0.48 

0.48 

0.47 

0.46 

0.45 

0.45 

0.45 

0.46 

0.47 

0.49 

0.54 

0.60 

b72 

0.51 

0.51 

0.51 

0.50 

0.49 

0.48 

0.48 

0.49 

0.51 

0.53 

0.57 

0.63 

b73 

0.54 

0.54 

0.54 

0.53 

0.52 

0-52 

0.51 

0.52 

0.54 

0.56 

0.60 

0.66 

b74 

0.57 

0.57 

0.57 

0.55 

0.55 

0.54 

0.54 

0.55 

0.56 

0.59 

0.63 

0.68 

b75 

0.59 

0.59 

0.59 

0.58 

0.57 

0.57 

0.57 

0.57 

0.59 

0.61 

0.65 

0.70 

b76 

0.61 

0.61 

0.61 

0.59 

0.59 

0.59 

0.58 

0.59 

0.60 

0.63 

0.67 

0.72 

b77 

0.63 

0.63 

0.62 

0.61 

0.60 

0.60 

0.60 

0-61 

0.62 

0.64 

0.68 

0.73 

b78 

0.64 

0.64 

0.63 

0.62 

0.62 

0.61 

0.61 

0.62 

0.63 

0.65 

0.69 

0.74 

b79 

0.65 

0.65 

0.64 

0.63 

0.63 

0.62 

0.62 

0.63 

0.64 

0.66 

0.70 

0.75 

b80 

0.66 

0.66 

0.65 

0.64 

0.63 

0.63 

0.63 

0.64 

0.65 

0.67 

0.71 

0.76 

b81 

0.66 

0.66 

0.65 

0.64 

0.64 

0.63 

0.63 

0.64 

0.65 

0.67 

0.71 

0.76 

b82 

0.66 

0.66 

0.65 

0.64 

0.64 

0.63 

0.63 

0.64 

0.65 

0.67 

0.71 

0.76 

b83 

0.66 

0.66 

0.65 

0.64 

0.64 

0.63 

0.63 

0.64 

0.65 

0.67 

0.71 

0.76 

b84 

0.66 

0.66 

0.65 

0.64 

0.63 

0.63 

0.63 

0.64 

0.65 

0.67 

0.71 

0.76 

b85 

0.65 

0.65 

0.64 

0.63 

0.63 

0.62 

0.62 

0.63 

0.64 

0.66 

0.70 

0.75 

b86 

0.64 

0.63 

0.63 

0.62 

0.61 

0.61 

0.61 

0.61 

0.63 

0.65 

0.69 

0.74 

b87 

0.63 

0.62 

0.62 

0.61 

0.60 

0.60 

0.60 

0.60 

0.62 

0.64 

0.68 

0.73 

b88 

0.62 

0.62 

0.61 

0.60 

0.59 

0.59 

0.59 

0.60 

0.61 

0.63 

0.67 

0.72 

b89 

0.61 

0.61 

0.60 

0.59 

0.58 

0.58 

0.58 

0.59 

0.60 

0.62 

0.66 

0.71 

b90 

0.60 

0.59 

0.59 

0.58 

0.57 

0.57 

0.57 

0.57 

0.59 

0.61 

0.65 

0.70 

Table  5  (cont). 


b61 

b62 

b63 

b64 

b65 

b66 

b67 

b68 

b69 

b70 

b71 

b72 

b61 

1.00 

b62 

0.99 

1.00 

b63 

0.11 

0.12 

1.00 

b64 

0.21 

0.23 

0.60 

1.00 

b65 

0.44 

0.46 

0.53 

0.67 

1.00 

b66 

-0.10 

-0.10 

0.00 

0.00 

0.00 

1.00 

b67 

-0.09 

-0.09 

-0.01 

0.00 

0.01 

0.33 

1.00 

b68 

-0.06 

-0.06 

0.00 

0.01 

0.03 

0.82 

0.35 

1.00 

b69 

0.69 

0.70 

0.22 

0.39 

0.74 

-0.01 

0.01 

0.04 

1.00 

b70 

0.72 

0.73 

0.20 

0.36 

0.72 

-0.01 

0.00 

0.04 

1.00 

1.00 

b71 

0.75 

0.75 

0.18 

0.34 

0.70 

-0.02 

0.00 

0.03 

0.99 

1.00 

1.00 

b72 

0.77 

0.78 

0.18 

0.34 

0.69 

-0.02 

-0.01 

0.03 

0.99 

1.00 

1.00 

1.00 

b73 

0.80 

0.80 

0.18 

0.34 

0.69 

-0.03 

-0.01 

0.02 

0.98 

0.99 

1.00 

1.00 

b74 

0.82 

0.82 

0.17 

0.33 

0.68 

-0.03 

-0.02 

0.02 

0.97 

0.98 

0.99 

1.00 

b75 

0.83 

0.83 

0.17 

0.33 

0.67 

-0.03 

-0.02 

0.02 

0.97 

0.98 

0.99 

0.99 

b76 

0.84 

0.84 

0.17 

0.33 

0.66 

-0.04 

-0.02 

0.01 

0.96 

0.97 

0.98 

0.99 

b77 

0.85 

0.85 

0.17 

0.32 

0.66 

-0.04 

-0.02 

0.01 

0.96 

0.97 

0.98 

0.99 

b78 

0.86 

0.86 

0.17 

0.32 

0.66 

-0.04 

-0.03 

0.01 

0.95 

0.96 

0.98 

0.98 

b79 

0.87 

0.87 

0.17 

0.32 

0.65 

-0.04 

-0.03 

0.01 

0.95 

0.96 

0.97 

0.98 

b80 

0.87 

0.87 

0.17 

0.32 

0.65 

-0.04 

-0.03 

0.01 

0.94 

0.96 

0.97 

0.98 

b81 

0.88 

0.87 

0.17 

0.32 

0.64 

-0.04 

-0.03 

0.01 

0.94 

0.95 

0.97 

0.98 

b82 

0.87 

0.87 

0.17 

0.32 

0.64 

-0.04 

-0.03 

0.01 

0.94 

0.95 

0.97 

0.98 

b83 

0.87 

0.87 

0.16 

0.31 

0.64 

-0.04 

-0.03 

0.01 

0.94 

0.95 

0.97 

0.98 

b84 

0.87 

0.87 

0.17 

0.32 

0.64 

-0.04 

-0.03 

0.01 

0.94 

0.95 

0.97 

0.98 

b85 

0.87 

0.86 

0.16 

0.31 

0.64 

-0.04 

-0.03 

0.01 

0.94 

0.95 

0.97 

0.98 

b86 

0.86 

0.86 

0.16 

0.32 

0.65 

-0.04 

-0.02 

0.01 

0.95 

0.96 

0.97 

0.98 

b87 

0.85 

0.85 

0.16 

0.32 

0.65 

-0.04 

-0.02 

0.01 

0.95 

0.96 

0.97 

0.98 

b88 

0.84 

0.84 

0.16 

0.32 

0.65 

-0.04 

-0.02 

0.01 

0.95 

0.96 

0.98 

0.98 

b89 

0.84 

0.84 

0.16 

0.32 

0.65 

-0.04 

-0.02 

0.02 

0.96 

0.97 

0.98 

0.99 

b90 

0.83 

0.83 

0.16 

0.32 

0.66 

-0.03 

-0.02 

0.02 

0.96 

0.97 

0.98 

0.99 

b91 

0.82 

0.82 

0.16 

0.31 

0.66 

-0.03 

-0.02 

0.02 

0.96 

0.97 

0.98 

0.99 

b92 

0.81 

0.81 

0.16 

0.31 

0.66 

-0.03 

-0.02 

0.02 

0.97 

0.98 

0.99 

0.99 

b93 

0.14 

0.14 

0.03 

0.06 

0.15 

-0.01 

0.22 

0.00 

0.21 

0.21 

0.21 

0.22 

b94 

0.78 

0.78 

0.15 

0.30 

0.64 

-0.03 

-0.01 

0.02 

0.94 

0.95 

0.96 

0.96 

b95 

0.26 

0.29 

0.35 

0.47 

0.68 

0.03 

0.03 

0.06 

0.71 

0.68 

0.66 

0.64 

b96 

-0.20 

-0.20 

-0.02 

-0.03 

-0.04 

0.41 

0,34 

0.51 

-0.07 

-0.08 

-0.09 

-0.09 

b97 

0.02 

0.02 

0.02 

0.06 

0.12 

0.00 

0.00 

0.01 

0.16 

0.16 

0.15 

0.15 

b98 

-0.15 

-0.15 

-0.02 

-0.01 

-0.01 

0.25 

0.34 

0.27 

-0.02 

-0.02 

-0,03 

-0.04 

b99 

0.11 

0.12 

0.06 

0.11 

0.23 

0.32 

0.01 

0.35 

0.31 

0.31 

0.30 

0.29 

b100 

0.47 

0.48 

0.17 

0.32 

0.68 

0.01 

0.03 

0.06 

0.95 

0.94 

0.93 

0.91 

b101 

0.49 

0.50 

0.17 

0.33 

0.69 

0.01 

0.03 

0.06 

0.95 

0.95 

0.94 

0.92 

b102 

0.52 

0.53 

0.17 

0.33 

0.69 

0.01 

0.03 

0.06 

0.96 

0.96 

0.94 

0.93 

b103 

0.54 

0.54 

0.17 

0.33 

0.68 

0.01 

0.02 

0.06 

0.96 

0.96 

0.95 

0.94 

b104 

0.55 

0.56 

0.17 

0.33 

0.69 

0.01 

0.02 

0.05 

0.97 

0.96 

0.95 

0.94 

b105 

0.56 

0.57 

0.17 

0.33 

0.68 

0.00 

0.02 

0.05 

0.97 

0.96 

0.96 

0.95 

b106 

0.58 

0.58 

0.17 

0.33 

0.68 

0.00 

0.02 

0.05 

0.97 

0.97 

0.96 

0.95 

b107 

0.59 

0.59 

0.17 

0.33 

0.69 

0.00 

0.02 

0.05 

0.97 

0.97 

0.97 

0.96 

b108 

0.60 

0.61 

0.17 

0.33 

0.69 

0.00 

0.02 

0.05 

0.98 

0.98 

0.97 

0.96 

b109 

0.62 

0.63 

0.17 

0.34 

0.69 

0.00 

0.01 

0.05 

0.98 

0.98 

0.98 

0.97 

b110 

0.63 

0.63 

0.17 

0.33 

0.69 

0.00 

0.01 

0.05 

0.98 

0.98 

0.98 

0.97 

b111 

0.63 

0.63 

0.17 

0.33 

0.69 

0.00 

0.01 

0.05 

0.98 

0.98 

0.98 

0.97 

b112 

0.62 

0.62 

0.17 

0.33 

0.69 

0.00 

0.02 

0.05 

0.98 

0.98 

0.97 

0.97 

b113 

0.60 

0.60 

0.17 

0.34 

0.69 

0.00 

0.02 

0.05 

0.97 

0.97 

0.97 

0.96 

b114 

0.58 

0.58 

0.17 

0.33 

0.69 

0.00 

0.02 

0.05 

0.97 

0.97 

0.96 

0.95 
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Table  5  (cont.). 


b85 

b86 

b87 

b88 

b89 

b90 

b92 

b93 

b94 

b95 

b96 

b85 

1.00 

b86 

1,00 

1.00 

b87 

1.00 

1.00 

1.00 

b88 

1.00 

1.00 

1.00 

1.00 

b89 

1.00 

1.00 

1.00 

1.00 

1.00 

b90 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

b91 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

b92 

0.99 

0.99 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

b93 

0,21 

0.21 

0.21 

0.21 

0.21 

0.21 

0.21 

0.21 

1.00 

b94 

0.96 

0.96 

0.97 

0.97 

0.97 

0.97 

0.97 

0.97 

0.21 

1.00 

b95 

0.54 

0.55 

0.56 

0.56 

0.57 

0.57 

0.58 

0.58 

0.14 

0.56 

1.00 

b96 

-0.12 

-0.11 

-0.11 

-0.11 

-0.11 

-0.11 

-0.10 

-0.10 

-0.03 

-0.10 

0.00 

1.00 

b97 

0.12 

0.12 

0.12 

0.13 

0.13 

0.13 

0.13 

0.14 

0.03 

0.13 

0.15 

0.14 

b98 

-0.06 

-0.06 

-0.06 

-0.06 

-0.05 

-0.05 

-0.05 

-0.05 

0.17 

-0.04 

0.03 

0.32 

b99 

0.25 

0.26 

0.26 

0.26 

0.26 

0.27 

0.27 

0.27 

0.06 

0.27 

0.26 

0.19 

b100 

0.82 

0.83 

0.84 

0.84 

0.85 

0.86 

0.87 

0.88 

0.20 

0.86 

0.73 

-0.03 

b101 

0.84 

0.84 

0.85 

0.86 

0.87 

0.87 

0.88 

0.89 

0.20 

0.87 

0.74 

-0.03 

b102 

0.85 

0.86 

0.87 

0.87 

0.88 

0.89 

0.89 

0.90 

0.21 

0.88 

0.74 

-0.04 

b103 

0.86 

0.87 

0.88 

0.88 

0.89 

0.90 

0.91 

0.91 

0.21 

0.89 

0.72 

-0.04 

b104 

0.87 

0.88 

0.89 

0.90 

0.90 

0.91 

0.91 

0.92 

0.21 

0.90 

0.72 

-0.04 

b105 

0.88 

0.89 

0.90 

0.90 

0.91 

0.92 

0.92 

0.93 

0.21 

0.90 

0.71 

-0.04 

b106 

0.89 

0.90 

0.91 

0.91 

0.92 

0.92 

0.93 

0.93 

0.21 

0.91 

0.71 

-0.05 

b107 

0.90 

0.91 

0.91 

0.92 

0.92 

0.93 

0.94 

0.94 

0.21 

0.92 

0.71 

-0.05 

b108 

0.91 

0.92 

0.92 

0.93 

0.93 

0.94 

0.94 

0.95 

0.21 

0.92 

0.70 

-0.05 

b109 

0.92 

0.93 

0.93 

0.94 

0.94 

0.95 

0.95 

0.95 

0.21 

0.93 

0.70 

-0.06 

b110 

0.92 

0.93 

0.93 

0.94 

0.94 

0.95 

0.95 

0.96 

0.21 

0.93 

0.70 

-0.06 

b111 

0.92 

0.93 

0.93 

0.94 

0.94 

0.95 

0.95 

0.95 

0.21 

0.93 

0.70 

-0.06 

b112 

0.91 

0.92 

0.93 

0.93 

0.94 

0.94 

0.95 

0.95 

0,21 

0.93 

0.70 

-0.05 

b113 

0.90 

0.91 

0.92 

0.92 

0.93 

0.93 

0.94 

0.94 

0.21 

0.92 

0.71 

-0.05 

b114 

0.89 

0.90 

0.91 

0.91 

0.92 

0.92 

0.93 

0.93 

0.21 

0.91 

0.71 

-0.05 

b115 

0.88 

0.89 

0.90 

0.90 

0.91 

0.91 

0.92 

0.92 

0.21 

0.90 

0.72 

-0.04 

b116 

0.00 

0.00 

0.00 

0.01 

0.01 

0.01 

0.24 

0.01 

0.05 

0.21 

b117 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

0.00 

0.02 

0.06 

0.15 

b118 

-0.08 

-0.07 

-0.07 

-0.07 

-0.07 

-0.06 

-0.06 

-0,06 

0,16 

-0.06 

0.01 

0.32 

b119 

-0.14 

-0.14 

-0.13 

-0.13 

-0.13 

-0.13 

-0.12 

-0.12 

0.10 

-0.11 

-0.01 

0.36 

b120 

-0.12 

-0.11 

-0.11 

-0.11 

-0.10 

-0.10 

-0.10 

0.12 

-0.09 

0.01 

0.31 

b121 

-0.18 

-0.17 

-0.17 

-0.17 

-0.17 

-0.16 

-0.16 

-0.16 

0.07 

-0.15 

-0.03 

0.41 

b122 

-0.14 

-0.14 

-0.13 

-0.13 

-0.13 

-0.12 

0.11 

-0.11 

-0.01 

0.33 

b123 

-0.17 

-0.17 

-0.16 

-0.16 

-0.16 

-0.16 

-0.16 

-0.16 

0.08 

-0.15 

-0.03 

0.32 

b124 

-0.20 

-0.19 

-0.19 

-0.19 

-0.19 

-0.18 

-0.18 

-0.18 

0.10 

-0.17 

-0.07 

0.38 

b125 

-0.24 

-0.23 

-0.23 

-0.23 

-0.23 

-0.23 

-0.23 

-0.23 

-0.03 

-0.23 

-0.12 

0.07 

b126 

0.44 

0.46 

0.46 

0.47 

0.47 

0.48 

0.48 

0.48 

0.11 

0.47 

0.56 

0.02 
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Table  6. 


Selected  latent  roots  of  the  correlation  matrix  for  the 
hyperspectral  image  data 


Latent 

Roots 

Latent  roots 

%  variance 

1 

76.8 

60.9 

2 

23.7 

18.8 

3 

4.55 

3.61 

4 

4.10 

3.25 

5 

1.83 

1.45 

6 

1.74 

1.38 

7 

1.31 

1.04 

8 

1.26 

1.00 

9 

1.01 

0.80 

10 

0.92 

0.73 

This  part  of  the  analysis  is  at  an  early  stage.  We  are  now  exploring  ways  of  improving  the  PCA  to 
obtain  more  refined  information  about  which  bands  are  likely  to  provide  relevant  spatial 
information.  The  correlations  between  the  original  wavebands  and  the  principal  components  have 
been  computed  and  those  bands  with  correlations  greater  than  0.9  on  the  first  5  components  only 
have  been  retained  for  a  further  PCA.  This  part  of  the  work  will  be  given  in  detail  in  the  next 
interim  report. 


Table  8.  Fitted  variogram  models  for  principal  components  (PC)  1  to  4. 


Variable 

Model 

Model  parameters 

^0 

Cl 

C2 

ai 

dl 

PC  1 

Double  spherical 

25.32 

19.02 

33.18 

3.716 

20.36 

PC  2 

Double  spherical 

8.318 

6.918 

7.536 

4.944 

14.26 

PC  3 

Double  spherical 

1.951 

1.445 

1.335 

4.683 

10.81 

PC  4 

Double  spherical 

1.940 

0.8951 

0.7950 

3.939 

13.67 

Note:  Co  is  the  nugget  variance  C]  and  C2  the  sills  of  the  autocorrelated  variance  and  a,  and  02  the  range 


of  spatial  dependence. 
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Appendix 


tn  c  ^  n  n  <  ,  r  o  ,  n  n  o 


AC  ****  PROGRAM  TO  COMPUTE  MOVING  AVERAGES  AND  VARIANCES  FOR  SQUARES 
OF  VARIOUS  SIZES 
^  ****  R  WEBSTER  ROTHAMSTED 
C  Latest  version  22  July  1999 

This  program  was  written  as  part  of  US  project 
C  and  may  be  handed  over  to  TEC . 

Program  reads  data  on  a  grid  with  X  and  Y  coordinates 
and  converts  them  to  an  array  for  the  selected  variate 
with  implied  coordinates. 

DIMENSION  ZK(190, 189) ,  ZA(30),  VM(190,189),  AM(190,189) 

****  ZK(  ,  )  will  hold  grid  of  data. 

character*72  TITLE (2) 

character*72  INFILE,  0P12,  INll,  FDAT 
data  MAXROW,  MAXCOL/190 , 189/ 
data  IN, INDAT,LP/10, 11, 12/ 

PRINT  *  ,  'WHAT  IS  THE  NAME  OF  THE  STEERING  FILE  ?' 

READ  (5, ' (A) ' )  INFILE 

OPEN  (INDAT, FILE=INFILE, STATUS= 'OLD' ) 
print  *  ,  'WHAT  IS  THE  NAME  OF  THE  DATA  FILE  ?  ' 

read  (5, ' (a) ' )  INll 
open  ( IN, f ile=INll , status= ' OLD ' ) 

PRINT  *, 'WHAT  DO  YOU  WANT  TO  CALL  THE  MAIN  OUTPUT  FILE  ?' 

READ  (5,  '  (A)  '  )  OP12 
OPEN  (LP, FILE=OP12 , STATUS= 'NEW' ) 

PRINT  * , ' WHAT  DO  YOU  WANT  TO  CALL  THE  SECOND  RESULTS  FILE  ? ' 
READ  (5,  '  (A)  '  )  0P12 
OPEN  (LF8 , FILE=0P12 , STATUS= ' NEW ' ) 

READ  (INDAT, 10)  TITLE 
WRITE  (LP,10)  TITLE 
WRITE  (LF8,12)  TITLE 
10  FORMAT  (A) 

NVAR  =  int ( CYNPUT ( INDAT ) +  0 . 1 ) 

NSEL  =  int (CYNPUT (INDAT) +0 . 1) 

MSIDE  =  int (CYNPUT (INDAT) +0 . 1) 

ZMIS  =  CYNPUT (INDAT) 

ILOG  =  int (CYNPUT (INDAT) +0.1) 
if  (ILOG.eq.l)  SHIFT=CYNPUT ( INDAT) 

****  NVAR  is  number  of  variates  in  file. 

NSEL  is  the  one  selected  for  analysis. 

MSIDE  is  the  side  of  the  square  within  which 
averages  are  computed. 

ZMIS  is  the  value  used  for  missing  or  blank. 

ILOG  =  1  to  transform  to  log  to  base  10 . 

SHIFT  is  a  value  to  be  added  to  data  to  shift  the  origin 
before  taking  logarithms. 

****  Set  data  grid  to  blank 

if  (ILOG.eq.l)  ZMIS=logl0 (ZMIS) 
do  20  1=1, MAXROW 
do  20  J=l,MAXCOL 

ZK(I, J) =ZMI S- 10000 . 0 
20  continue 

read  (INDAT, 10)  FDAT 
^  ****  Read  the  data. 

35  NC=0 
NROW=0 
NCOL=0 

36  NC=NC+1 

read  (IN, FDAT, end=45)  ICOL,  IROW,  (ZA(J),  J=1,NVAR) 
if  (ICOL.gt .MAXCOL)  then 
write  (LP,38)  ICOL 


stop 

endif 

if  (IROW.gt.MAXROW)  then 
write  (LP,39)  IROW 
stop 
endif 

38  format  (/lOx, 'ICOL  exceeds  array  bound'//) 

39  format  (/lOx, 'IROW  exceeds  array  bound'//) 
if  (NROW.lt. IROW)  NROW=IROW 

if  (NCOL.lt. ICOL)  NCOL=ICOL 

ZL=ZA(NSEL) 

if  (ILOG.eq.l)  then 

if  (ZL.gt.0.01)  then 
ZL=loglO (ZL+SHIFT) 
else 

ZL=ZMIS-100000 . 0 
endif 
endif 

ZKdROW,  ICOL)  =ZL 
goto  36 
45  continue 
NC=NC-1 

write  (LP,47)  NC,  NROW,  NCOL 
47  format (//  lOx,  'Number  of  data  ',ilO/ 

1  lOx,  'Number  of  rows  ',ilO/ 

2  lOx,  'Number  of  columns  ',ilO/) 

if  (ILOG.eq.l)  write  (LP,  51)  SHIFT 

51  format  (/lOx, 'DATA  TRANSFORMED  TO  LOG  TO  BASE  10'/ 

1  lOx, 'SHIFT  ' ,F10.3/) 

ZMAX=-99999999 
ZMIN=999999999 
NN=NC 
ZBAR=0 . 0 
SSQ=0 . 0 
COUNT=0.0 
do  54  1=1, NROW 
do  53  J=l,NCOL 
ZZ=ZK(I, J) 

if  (ZZ.le.ZMIS)  goto  53 

if  (ZMAX.lt.ZZ)  ZMAX=ZZ 

if  (ZMIN.gt.ZZ)  ZMIN=ZZ 

DIF=ZZ-ZBAR 

COUNT=COUNT+l . 0 

ZBAR=ZBAR+DIF/COUNT 

SSQ=SSQ+(1.0-1.0/COUNT) *DIF*DIF 

53  continue 

54  CONTINUE 
A3  =  0 . 0 

do  57  1=1, NROW 
do  56  J=l,NCOL 
ZZ=ZK(I, J) 

if  (ZZ.le.ZMIS)  goto  56 
A3=A3+ (ZZ-ZBAR) **3 

56  continue 

57  CONTINUE 
A2=SSQ/ COUNT 

A3 = (A3 /COUNT) / (A2*sqrt (A2) ) 

VAR=SSQ/ (COUNT- 1.0) 

STD=sqrt (VAR) 

write  (LP,  58)  COUNT,  ZMIN,  ZMAX,  ZBAR,  VAR,  STD,  A3 

58  format  (//,10X,  '  Count  ',fl0.l/ 

1  lOx,  '  Minimum  ',fl0.4/ 

1  lOx,  '  Maximum  ',fl0.4/ 

1  lOx,  '  Mean  ',fl0.4/ 


2  lOx,  '  Variance  ',fl2.6/ 

3  lOx,  '  Standard  deviation  ',fl0.4/ 

4  lOx,  '  Skewness  ',fl0.4/) 

****  Enter  means  and  variances  into  arrays  VM (  ,  )  and  AM(  ,  ) . 
do  122  I=l,MAXROW 
do  122  J=l,MAXCOL 
AM(I, J) =ZBAR 
VM(I, J) =VAR 

122  continue 

****  Now  start  computing  moving  averages  and  moving  variances. 

Data  are  in  array  ZK(  ,  )  with  NCOL  columns  and  NROW  rows. 
M=MSIDE 
NRM=NR0W-M-1 
NCM=NCOL-M-l 
do  200  KR=1,NRM 
do  200  KC=1,NCM 

****  Initialize  counters 
COUNT=0 . 0 
SUM=0.0 
SSQ=0.0 

****  Loop  through  each  square  in  turn 
KRS=KR 
KRE=KR+M-1 
KCS=KC 
KCE=KC+M-1 
do  145  I=KRS,KRE 

do  145  J=KCS,KCE 
ZZ=ZK(I, J) 

if  (ZZ.lt.ZMIS)  goto  145 
COUNT=COUNT+l . 0 
DIF=ZZ-SUM 
SUM=SUM+DIF/ COUNT 
SSQ=SSQ+(1.0-1.0/COUNT) *DIF*DIF 

145  continue 

SSQ=SSQ/ (COUNT-1. 0) 

SDV=sqrt (SSQ) 

KOR=KRS+ (M+1) /2 
KOC=KCS+ (M+1) /2 
VM(KOR,KOC) =SSQ 
AM(KOR,KOC) =SUM 

200  continue 

write  (LP,233) 
do  224  1=1, NROW 

write  (LP,  230)  (AM(I,J),  J=l,NCOL) 

224  continue 
write  (LP, 238) 
do  225  1=1, NROW 

write  (LP,  231)  (VM(I,J),  J=l,NCOL) 

225  continue 

226  format  (//) 

230  format  (2x,f9.4) 

231  format  (2x,fl0.5) 

233  format  (//lOx,  'MOVING  AVERAGES  '/) 

238  format  (//lOx,  'MOVING  VARIANCES  '/) 
stop 
end 

FUNCTION  CYNPUT(IN) 

****  READS  A  REAL  NUMBER  FROM  AN  80 -BYTE  RECORD  IN  FREE  FORMAT 
DIMENSION  K(80) ,NUM(10) 

DATA  NUM/lHO , IHl , 1H2 , 1H3 , 1H4 , 1H5 , 1H6 , 1H7 , 1H8 , 1H9/ 

DATA  INOLD,N, IFL,NPLUS,MINUS,NDOT/0, 81, 0, 1H+, IH- , IH. / 
CYNPUT=-0 .0 


IF(INOLD.EQ. IN.AND.  N.LE.80)  GOTO  20 
5  IF(IFL.NE.O) -RETURN 
INOLD=IN 

READ(IN,10)  (K(I) , 1=1, 80) 

10  FORMAT ( 8 OAl) 

15  N=1 

20  IF(N.GT.80)  GOTO  35 
DO  30  I=N,80 
II=K(I) 

DO  25  J=l,10 

IF (II .EQ.NUM(J) )  GOTO  40 
25  CONTINUE 

IF (II .EQ. MINUS)  GOTO  40 
IF(II.EQ.NDOT)  GOTO  40 
IF(II .EQ.NPLUS)  GOTO  40 
30  CONTINUE 
35  GOTO  5 
40  SIGN=1.0 

IFdI.EQ. MINUS)  SIGN=-1.0 
IF ( II. EQ. MINUS  .OR.  II. EQ.NPLUS)  1=1+1 
IF(I.GT.80)  GOTO  60 
DO  55  N=I,80 
NN=K(N) 

IF(NN.EQ.NDOT)GOTO  70 
DO  45  J=l,10 

KK=J-1 

IF(NN.EQ.NUM(J) )  GOTO  50 
45  CONTINUE 
GOTO  65 

50  CYNPUT=10 . 0*CYNPUT+KK 
55  CONTINUE 
60  N=82 

65  CYNPUT=SIGN*CYNPUT 
RETURN 
70  I=N+1 

TENS=1.0 

IF(I.GT.80)  GOTO  90 
DO  85  N=I, 80 
NN=K(N) 

DO  75  J=l,10 

KK=J-1 

IF(NN.EQ.NUM(J) )  GOTO  80 
75  CONTINUE 
GOTO  65 

80  TENS=TENS*0 . 1 

CYNPUT=CYNPUT+TENS*KK 
85  CONTINUE 
90  N=82 
GOTO  65 
END 


♦c  ****  PROGRAM  TO  COMPUTE  VARIOGRAMS  FOR  SQUARES  OF  VARIOUS  SIZES 
( 

u  ****  R  WEBSTER  ROTHAMSTED 
C  Laresc  version  22  July  1999 

( 

(  This  program  was  written  as  part  of  US  project 

C  and  may  be  handed  over  co  TEC. 

(  Program  reads  data  on  a  grid  with  X  and  Y  coordinates 

and  converts  them  to  an  array  for  the  selected  variate 
C-  with  implied  coordinates. 

It  tiles  the  grid  into  non -over lapping  squares  of 
C  given  side.  Any  points  to  the  bottom  or  right 

(  of  the  grid  left  over  play  no  role. 

i 

DIMENSION  ZK(190,189),  GRID(30,30),  ZA(30) 
real  WT(30),  GAM(30),  WIAG(30) 

(  ****  2K(  ,  )  will  hold  grid  of  data, 

character* 72  TITLE (2) 

character* 72  INFILE,  OP12 ,  INll,  OPGAM,  FDAT 
data  MAXROW,  MAXCOL/190 , 189/ 
data  IN, INDAT, LP, LF8/10 , 11 , 12 , 13/ 

PRINT  *  ,  'WHAT  IS  THE  NAME  OF  THE  STEERING  FILE  ?' 

READ  (.‘5,'(A)')  INFILE 

OPEN  (INDAT, FILE=INFILE,STATUS=' OLD' ) 

print  *  ,  'WHAT  IS  THE  NAME  OF  THE  DATA  FILE  ?  ' 

read  (5, ' (a) ' )  INll 

open  (IN, file=INll, status= 'OLD' ) 

PRINT  * , ' WHAT  DO  YOU  WANT  TO  CALL  THE  MAIN  OUTPUT  FILE  ? ' 

READ  (5,  '  (A)  '  )  OP12 

OPEN  (LP, FILE=OP12 , STATUS= 'NEW ' ) 

PRINT  *, 'WHAT  DO  YOU  WANT  TO  CALL  THE  FILE  OF  VARIOGRAMS?' 

READ  (5,  '  (A)  '  )  OPGAM 

OPEN  (LF8 , FILE=OPGAM, STATUS= 'NEW ' ) 

( 

READ  (INDAT, 10)  TITLE 
WRITE  (LP,10)  TITLE 
C  WRITE  (LF8,12)  TITLE 

10  FORMAT  (A) 

NVAR  =  int (CYNPUT ( INDAT) +0.1) 

NSEL  =  int (CYNPUT (INDAT) +0. 1) 

MSIDE  =  int (CYNPUT (INDAT) +0.1) 

MAXLAG=  int (CYNPUT ( INDAT) +0.1) 

ZMIS  =  CYNPUT (INDAT) 

ILOG  =  int (CYNPUT (INDAT) +0 . 1) 
if  (ILOG.eq.l)  SHI FT = C YNPUT ( INDAT) 

C  ****  NVAR  is  number  of  variates  in  file. 

(  NSEL  is  the  one  selected  for  analysis. 

(  MSIDE  is  the  side  of  the  square  within  which 

C  averages  are  computed. 

c  MAXLAG  is  the  maximum  lag  distance  of  variograms 

<  ZMIS  is  the  value  used  for  missing  or  blank. 

C  ILOG  =  1  to  transform  to  log  to  base  10. 

C  SHIFT  is  a  value  to  be  added  to  data  to  shift  the  origin 

!  before  taking  logarithms. 

^  ****  Set  data  grid  to  blank 
C  if  (ILOG.eq.l)  ZMIS=logl0 (ZMIS) 

do  20  1=1, MAXROW 
do  20  J=l,MAXCOL 

ZK(I, J) =ZMIS-10000 . 0 
20  continue 


read  (INDAT, 10)  FDAT 


****  Read  the  data. 

35  NC=0 
NROW=0 
NCOL=0 

36  NC=NC+1 

read  (IN, FDAT, end=45)  ICOL,  IROW,  (ZA(J) ,  J=1,NVAR) 
if  (ICOL.gt .MAXCOL)  then 
write  (LP,38)  ICOL 
stop 
endif 

if  (IROW.gt.MAXROW)  then 
write  (LP,39)  IROW 
stop 
endif 

38  format  (/lOx, 'ICOL  exceeds  array  bound'//) 

39  format  (/lOx, 'IROW  exceeds  array  bound'//) 
if  (NROW.lt .IROW)  NROW=IROW 

if  (NCOL.lt. ICOL)  NCOL=ICOL 

ZL=ZA(NSEL) 

if  (ILOG.eq.l)  then 

if  (ZL.gt.0.01)  then 
ZL=loglO (ZL+SHIFT) 
else 

ZL=ZMIS-100000 . 0 
endif 
endif 

ZKdROW,  ICOL)  =ZL 
goto  36 
45  continue 
NC=NC-1 

write  (LP,47)  NC,  NROW,  NCOL 
47  format (//  lOx,  'Number  of  data  ',ilO/ 

1  lOx,  'Number  of  rows  ',ilO/ 

2  lOx,  'Number  of  columns  ',ilO/) 

if  (ILOG.eq.l)  write  (LP,  51)  SHIFT 

51  format  (/lOx, 'DATA  TRANSFORMED  TO  LOG  TO  BASE  10'/ 

1  lOx, 'SHIFT  ' ,F10.3/) 

ZMAX=:-99999999 
ZMIN=999999999 
NN=NC 
ZBAR=0.0 
SSQ=0. 0 
COUNT=0.0 
do  54  1=1, NROW 
do  53  J=l,NCOL 
ZZ=ZK(I, J) 

if  (ZZ.le.ZMIS)  goto  53 

if  (ZMAX.lt.ZZ)  ZMAX=ZZ 
if  (ZMIN.gt.ZZ)  ZMIN=ZZ 
DIF=ZZ-ZBAR 
COUNT=COUNT+l . 0 
ZBAR=ZBAR+DIF/ COUNT 
SSQ=SSQ+(1.0-1.0/COUNT) *DIF*DIF 

53  continue 

54  CONTINUE 
A3  =  0 . 0 

do  57  1=1, NROW 
do  56  J=l,NCOL 
ZZ=ZK(I,J) 

if  (ZZ.le.ZMIS)  goto  56 
A3=A3+ (ZZ-ZBAR) **3 

56  continue 

57  CONTINUE 
A2=SSQ/COUNT 


A3= (A3 /COUNT) / (A2*sqrt (A2) ) 
VAR=SSQ/ (COUNT-1 . 0) 

STD=sqrt (VAR) 

write  (LP,  58)  COUNT,  ZMIN, 
58  format  (//,10X, 

1  lOx, 

1  lOx, 

1  lOx, 

2  lOx, 

3  lOx, 

4  lOx, 


Count 

Minimum 

Maximum 

Mean 

Variance 

Standard  deviation 
Skewness 


ZMAX,  ZBAR,  VAR,  STD, 
'  ,fl0.l/ 

'  , flO.4/ 

' ,fl0.4/ 

'  ,fl0.4/ 

' , fl2 . 6/ 

' ,fl0.4/ 

' ,fl0.4/) 


A3 


****  Compute  starting  in  top  left  corner  of  grid. 
NTILER=int (NROW/MSIDE) 

NTILEC=int (NCOL/MSIDE) 


do  300  IR=1,NTILER 
IRS= (IR-1) *MSIDE+1 
IRE=IR*MSIDE 
do  300  IC=1,NTILEC 
ICS= (IC-1) *MSIDE+1 
ICE=IC*MSIDE 
11  =  0 

do  210  I=IRS,IRE 
JJ=0 
11=11+1 

do  210  J=ICS,ICE 
JJ=JJ+1 

GRIDdI,  JJ)  =ZK(I,  J) 

210  continue 

****  Data  are  now  transferred  into  array  GRID(  ,  )  covering 
a  small  square  of  side  MSIDE. 

Initialize  accumulators, 
do  220  I=1,MAXLAG 
WLAG(I) =0 . 0 
GAM  (I)  =0 . 0 
WT(I) =0.0 
SUM=0 . 0 


SSQ=0 . 0 
COUNT=0 . 0 


220  continue 

do  225  1=1, MSIDE 
do  225  J=l, MSIDE 
ZZ=GRID(I, J) 

if  (ZZ.lt.ZMIS)  goto  225 
COUNT=COUNT+l . 0 
DIF=ZZ-SUM 
SUM=SUM+DIF/ COUNT 
SSQ=SSQ+(1.0-1.0/COUNT) *DIF*DIF 
225  continue 

SSQ=SSQ/ (COUNT- 1.0) 

SDV=sqrt (SSQ) 

write  (LP,  227)  IRS,  ICS,  SUM,  SSQ,  SDV 
227  format  (//5X,  'CORDINATES  ',  16,16/ 

1  5X,  'MEAN  ',  F12.5/ 

2  5X,  'VARIANCE  ',  F12.5/ 

3  5X,  'ST.  DEVIATION  ',  F12.5/) 

write  (LP,230) 

230  format  (/2x,  'LAG  SEMI VARIANCE  COUNT'/) 

write  (LF8,234)  IRS, ICS, IRS, ICS, SSQ 
234  format  (lx, ' COORDINATES ', 2i6 , '  XY ' / ' : ' /2x, 2i8 , F12 . 6/ ' : ' ) 
do  255  1=1, MSIDE 
do  255  J=l, MSIDE 
Z1=GRID(I,J) 


if  (Zl.lt.ZMIS)  goto  255 
do  245 •K=1,MSIDE 
do  245  L=1,MSIDE 
Z2=GRID (K,L) 

if  (Z2.1t.ZMIS)  goto  245 
Xl  =  float  (I) 

X2=float (K) 

Yl=float (J) 

Y2=float (L) 

D=sqrt ( (X1-X2) **2+ (Y1-Y2) **2) 

LAG=int (D) +1 
WLAG (LAG) =WLAG (LAG) +D 
GAM(LAG) =GAM(LAG) + (Z1-Z2) **2 
WT (LAG) =WT (LAG) +1 . 0 
245  continue 

255  continue 
ANGLE=0 . 0 
do  270  I=1,MAXLAG 

WLAG(I) =WLAG(I) /WT(I) 

GAM(I)=0.5*GAM(I) /WT(I) 

write  (LP,275)  WLAG (I) , GAM (I) , WT (I) 

write  (LF8,275)  WLAG(I) ,GAM(I) ,WT (I) 

270  continue 

275  format  (2x, f 7 . 2 , f 12 . 6 , f 10 . 1 ) 
write  (LF8,280) 

280  format  ( ' : ' ) 

300  continue 
stop 
end 

FUNCTION  CYNPUT(IN) 

****  READS  A  REAL  NUMBER  FROM  AN  80 -BYTE  RECORD  IN  FREE  FORMAT 
DIMENSION  K(80)  ,NUM(10) 

DATA  NUM/ IHO , IHl , 1H2 , 1H3 , 1H4 , 1H5 , 1H6 , 1H7 , 1H8 , 1H9/ 

DATA  INOLD,N, IFL, NPLUS , MINUS , NDOT/0 , 81 , 0 , 1H+ , 1H-, IH./ 
CYNPUT=-0 . 0 

IFdNOLD.EQ.IN.AND.  N.LE.80)  GOTO  20 
5  IF(IFL.NE.O)  RETURN 
INOLD=IN 

READ(IN,10)  (K(I) , 1=1, 80) 

10  FORMAT ( 8 OAl) 

15  N=1 

20  IF(N.GT.80)  GOTO  35 
DO  30  I=N,80 
II=K(I) 

DO  25  J=l,10 

IF(II .EQ.NUM(J) )  GOTO  40 
2  5  CONTINUE 

IF(II .EQ. MINUS)  GOTO  40 
IF(II .EQ.NDOT)  GOTO  40 
IF(II .EQ.NPLUS)  GOTO  40 
30  CONTINUE 
35  GOTO  5 
40  SIGN=1.0 

IF(II .EQ. MINUS)  SIGN=-1.0 
IF(II .EQ. MINUS  .OR.  II. EQ.NPLUS)  1=1+1 
IF(I .GT.80)  GOTO  60 
DO  55  N=I , 80 
NN=K(N) 

IF (NN. EQ.NDOT) GOTO  70 
DO  45  J=l,10 
KK=J-1 

IF(NN.EQ.NUM(J) )  GOTO  50 
45  CONTINUE 


GOTO  65 

5  0  CYNPUT= 10.0*  CYNPUT+KK 

55  CONTINUE 
60  N=82 

65  CYNPUT=SIGN*CYNPUT 
RETURN 
70  I=N+1 
TENS=1 . 0 

IF(I.GT.80)  GOTO  90 
DO  85  N=I,80 
NN=K{N) 

DO  75  J=l,10 
KK=J-1 

IF(NN.EQ.NUM(J) )  GOTO  80 
75  CONTINUE 
GOTO  65 

80  TENS=TENS*0 . 1 

CYNPUT=CYNPUT+TENS*KK 
85  CONTINUE 
90  N=82 
GOTO  65 
END 


